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Up to date, the life time of experimentally demonstrated entangled states has been limited, due to 
their fragility under decoherence and dissipation. Therefore, they are created under strict isolation 
conditions. In contrast, new approaches harness the coupling of the system to the environment, 
which drives the system into the desired state. Following these ideas, we present a robust method 
for generating steady state entanglement between two distant atomic ensembles. The proposed 
scheme relies on the interaction of the two atomic systems with the common vacuum modes of 
the electromagnetic field which act as an engineered environment. We develop the theoretical 
framework for two level systems including dipole-dipole interactions and complement these results 
by considering the implementation in multi-level ground states. 

PACS numbers: 03.67.Mn, 32.80.Qk 
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I. INTRODUCTION 

Entanglement is the most prominent and distinctive 
feature of quantum mechanics. As it plays a central role 
in fundamental tests of Quantum Mechanics and in ap- 
plications in the field of quantum information science, 
entanglement has been generated and studied in various 
systems. However, experimental generation of entangled 
states and their use in quantum information and commu- 
nication protocols is hampered by their fragility under 
decoherence. The life time of entangled quantum states 
is typically very short. The pursuit of the generation of 
persistent entanglement is not only of fundamental in- 
terest in view of the investigation of entangled quantum 
states at long time-scales but also vital for many applica- 
tions, such as quantum repeaters [THTT]. On account of 
this problem, quantum systems are usually strictly iso- 
lated in the endeavor to avoid their interaction with the 
environment. 

In contrast, we adopt an ostensibly counter-intuitive 
approach using dissipation [T2ti46] . Here, the interaction 
of the system with the environment is employed such that 
dissipation drives the system into the desired state. More 
specifically, we propose and analyze a scheme for the gen- 
eration of long-lived entanglement between two distant, 
mesoscopic ensembles (see also |17j). Both atomic sam- 
ples are placed in magnetic fields and interact with an 
environment consisting of the vacuum modes of the elec- 
tromagnetic field. A laser field mediates the coupling of 
the atomic system to the environment. The interaction 
of the system and the bath can be controlled via laser- 
and magnetic fields, which allow one to engineer the cou- 
pling in such a way that the unique steady state of the 
dissipative evolution is an entangled state. 

This dissipative approach has several remarkable ad- 
vantages. For example, the scheme performs well start- 
ing from an arbitrary initial state. This feature renders 
the initialization of the system in a pure state unneces- 
sary. Most importantly, the evolution is robust against 
moderate external noise. Entanglement is obtained in a 



steady state. This auspicious property is very promis- 
ing in view of the quest for viable, extremely long-lived 
entanglement. 

We develop a scheme for two-level systems and show 
that steady-state entanglement can be generated in the 
presence of undesired transitions and fluctuating mag- 
netic fields. We also include external pump fields and 
find that surprisingly, incoherent pumping can be ben- 
eficial in a certain parameter range. These central re- 
sults are supplemented by a short and very general study 
of the implementation in atoms with muli-level ground 
states. Since additional dynamics in atoms with multi- 
level structure lead to particle losses, a quasi-steady state 
is produced. Remarkably, incoherent pumping enables 
the creation of entanglement in a true steady state. We 
investigate the conditions for the generation of long-lived 
entanglement in mesoscopic multi-level systems, using a 
simplified model and consider the realization of the pro- 
posed method in ^^'^Cs vapors as used in P51 - [5U| as spe- 
cific example. 

The paper is organized as follows. The main idea is 
introduced in Sec. IHJ which also contains a summary of 
the central results. In Sec. IIIII we derive the full mas- 
ter equation for two-level atoms and calculate the evo- 
lution of entanglement. Prospects for generating steady- 
state entanglement in multi-level systems are discussed 
in Sec. |IV| Sec. \V\ summarizes and concludes the paper. 



II. MAIN IDEA AND CENTRAL RESULTS 

In the following, we explain the basic idea for gener- 
ating purely dissipatively driven entanglement in atomic 
ensembles. We introduce a realistic description including 
noise effects and discuss the prospects for realizing the 
proposed scheme experimentally. 

We start by explaining the underlaying concept for two 
bosonic modes with annihilation operators a and b. The 
entangled target state under consideration is a two mode 
squeezed (TMS) state |^tms)i which is characterized in 




FIG. 1: (Color online) Setup for creating steady-state entan- 
glement between two atomic ensembles separated by a dis- 
tance R. Both ensembles are placed in magnetic fields B, 
which are oriented along x. A strong y-polarization laser 
beam (shown in green color) propagates along z and couples 
the atomic system to the environment which consists of the 
vacuum modes of the copropagating electromagnetic field in x 
polarization (depicted in red). The interaction between atoms 
and light modes is illustrated in Fig. [2] 
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where the nonlocal annihilation operators A and B 
are given by 

A = ^1 a + V h\ 
B = /i 6 + i^ a ' . 



(1) 



This equation completely characterizes a particular 
squeezed state with squeezing parameter r where ^ — 
cosh(r) and v = sinh(r). This state can be prepared by 
means of a dissipative evolution governed by the master 
equation 

dtp{t) = K^ [Ap{t)A^ - A^Ap{t)/2 - p{t)A^A/2^ 

+ Kg {Bp{t)B^ - B'^Bp{t)/2 - p(t)B'fB/2] , (2) 

where the coefficients k^ and kj^ are positive. This time 
evolution drives the system into the unique steady state 
Poo = |*TMs)(^TMs|, for i ^ oo (see App.JAJ). 
Starting from this result, we explain how this concept 
can be applied for creating entanglement between two 
macroscopic atomic ensembles and how a dissipative evo- 
lution of type (pi) can be realized. We consider two atomic 
ensembles placed in magnetic fields as shown in Fig. IT] 
Atoms are assumed to possess a two-level ground state 
with internal states |t) and \l) (see Fig. [2]). We aim at 
creating an entangled state, where the collective spins of 
the two samples are correlated in y and in z direction. 
The amount of entanglement generated can be measured 



by the quantity 



var(Jyj + Jyji) -I- var(Jz,i - Jz.ii) 
\{J.,i)\ + \{J.Ji)\ 



(3) 
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For separable states ^ > 1. Jx,i = J2i=ifx i i^ ^^^ '^^^~ 
lective spin in x direction in the first ensemble, where 
Ni is the number of atoms and jl j is the x-spin compo- 
nent of the zth atom [S3] . Analogous definitions hold for 
for spins in y and z direction and for the collective spin 
of the second ensemble. In order to prepare this target 
state, we use a dissipative evolution governed by a mas- 
ter equation of type ([2| , where the nonlocal operators A 
and B are replaced by 
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and ^7/77 are collective spin operators 54i, with 

^Ef=ili).(t|such 
that Jy = i(J+ + J-) and J^ == |(J+ - J-). Nor- 
malization of the operators [A^ A^ = [B, B^ = 1 re- 
quires p^ — v^ = 1. The light-matter interaction shown 
in Fig. [2] gives rise to the desired master equation. More 
specifically, after adiabatic elimination of excited states, 
the effective ground state Hamiltonian is of the form 

H « /a.„ d\. (Aal + A^a^) + J^^^^ dk (Saj, + BU 
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where a^. is the creation operator for a photon with wave 
vector k. The first and second integral cover narrow 
bandwiths Awj^ and Aujus centered around the lower 



and upper sideband respectively (see Sec. Ill A). The 



modes of the light field are treated as bath and are 
therefore traced out. Using the Born-Markov approxima- 
tion, one obtains a master equation of standard Lindblad 
form (compare Eq. (|2|). Collective Lamb shifts can be 
shown to be negligible in the setting considered here (see 
App. B2). In the limit t — >■ 00, this evolution drives the 



system into an entangled steady state. In the absence of 
noise and for A^ ^ 1, 
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Next, we include additional processes such as thermal 
motion, undesired atomic transitions and fluctuating 
magnetic fields as well as resonant pump fields. Details 
can be found in Sec. |III C[ For large particle numbers 
Nj = Nil = A^ > 1 and t -> 00, we find 
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where ^2,00 is the steady state value of the atomic 
polarization P2 = ^ Eti (It).Jtl - ID^UI), T is the 
single particle decay rate and F is the dephasing rate 
associated with noise effects, d is the optical depth 
of one atomic ensemble. As shown in Sec. |IIID[ the 
application of resonant pump fields can be beneficial 



even though noise is added by doing so. Note that 
for d — > oo, ^oo — >■ C^°^'- For a large optical depth, 
the entangling dynamics is significantly enhanced by 
collective effects. In contrast, noise processes are single 
particle effects and therefore not amplified by a factor d. 
Eq. (51) shows that for strong coupling between atoms 
and light, entanglement can be generated in a steady 
state. This is the main result of this article. 

Intuitively, entanglement is created by virtue of in- 
terference of different processes in the first and second 
ensemble. As illustrated in Fig.[2j the interaction of light 
and atoms in the first ensemble is chosen such that the 
emission of a photon in the upper sideband corresponds 
to a spin flip |t) —>■ \l) . Similarly, the emission of a pho- 
ton in the upper sideband involves a spin flip \\,) — >■ jf) in 
the second ensemble. Due to collective effects [H], light 
is emitted in forward direction with high probability, 
for hot samples with high optical depth. As spin flips 
in either ensemble lead to emission of light into the 
same spatial mode, both processes are indistinguishable 
if a photon is detected. (An analogous argument holds 
for photons in the lower sideband.) In this respect, 
the setup resembles quantum repeater schemes [IHin], 
where collective excitations in two atomic ensembles 
are converted to photons, which subsequently interfere 
at a 50/50 beamspitter such that entanglement can be 
created conditioned on the detection of a photon in 
one of the two output ports of the beamsplitter |55j . 
Here, no beamsplitter is needed, since both ensembles 
emit into the same spatial mode. The most important 
difference, however, lies in the fact that our scheme is 
not conditioned on a specific measurement outcome. It 
works deterministically and does not require a detection 
of the emitted photon, as the measurement is performed 
continuously by the environment. 
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FIG. 2: (Color online) Atomic level schemes, a) A magnetic 
field which is oriented along x (see Fig. [l| causes a Zeeman 
splitting Q of atomic ground-states levels 1 1) and 1 1) and 
defines the quantization axis. A strong y polarized coherent 
field with detuning A drives transitions |t) — >■ lej) and \\) — >■ 
|e-f^). Coupling to the vacuum modes of the electromagnetic 
field gives rise to transitions le^.) — >■ \\) and |e-f-) — >■ |t). b) 
A static electric field is applied to the second ensemble such 
that the energy difference between ground and excited states 
is enhanced by 2 A. 



These fields add noise to the system and limit therefore 
the amount of entanglement that can be generated. 
However, for samples with high optical depth, incoher- 
ent pumping can render the creation of steady state 
entanglement in atoms with multi-level ground states 
possible. This is illustrated in Sec. |IVB| by considering 
^^•^Cs vapors at room temperature and experimental 
parameters close to the values published in [33 EOl ES] • 
We take the most fundamental limitations imposed by 
undesired radiative processes into account and estimate 
that steady state entanglement with ^oo = 0.9 should be 
attainable for a moderate optical depth d = 30 in the 
absence of other (implementation-dependent) sources of 



The ideas put forward in this work are devised 
and elaborated for two- levels systems, but the proposed 
scheme can also be realized using atoms with multi-level 
ground states. The scheme put forward in this work has 
been realized recently [SS] using alkai atoms. In this 
experiment, purely dissipatively driven entanglement 
between two macroscopic atomic ensembles at room 
temperature has been demonstrated yielding an order of 
magnitude improvement in the entanglement life time 
compared to previous experiments, where entanglement 
has been generated in this system using standard 
methods. In a multi-level setting, the population in the 
two-level subsystem is continually reduced due to unde- 
sired transitions to other ground state levels. We take 
this and other features of the multilevel structure into 
account and find that the continuous reduction of the 
collective spin leads to the production of a quasi steady 
state: the steady state with respect to the relevant 
two-level subsystem is superposed by slow additional 
dynamics due to the multi-level structure. This can be 
counteracted by the application of external pump fields. 



III. CREATION OF STEADY STATE 

ENTANGLEMENT IN A TWO-LEVEL SYSTEM 

As outlined above, light modes act as environment and 
the interaction between the system and the reservoir is 
controlled by means of laser- and magnetic fields. In 
Sec. Ill A[ we explain the interaction between atoms and 
light in more detail. The master equation governing the 
dissipative evolution of the reduced density matrix of the 
atomic system p is given by 



dtP^ ^cntP + ^n 



oP, 



where £ont and >Cnoise are Lindblad operators. Desired 
interactions give rise to the entangling dynamics repre- 
sented by ContP- The second term CnoiscP summarizes all 
undesired effects. 

Below, CentP and CnoiscP are determined. To this end, 
the master equation corresponding to the light-matter in- 
teraction in Figs.[l]and[2]is derived in Sec. HI B including 



undesired radiative processes. Excited states are adiabat- 
ically eliminated such that an effective master equation 
for atomic ground states 1 1) and 1 1) is obtained. In 
Sec. |III C[ thermal motion of atoms is taken into account 
and additional effects due to pump fields and noise pro- 
cesses are included. Based on these results, the amount 
of entanglement that can be produced is calculated in 
SeciniDl 



A. Light— matter interaction 

In this subsection, we describe the setup for creation 
of entanglement between two atomic ensembles and 
explain how light and matter interact. 

We consider the setup shown in Fig. [T] A strong 
y-polarized laser beam propagates along z and passes 
two atomic ensembles in a homogeneous magnetic field, 
which defines the quantization axis and is oriented along 
X. Each atomic ensemble consists of a large number N 
of hydrogen-like atoms with an internal level structure 
as depicted in Fig. [2] The laser field is assumed to cover 
a very narrow bandwidth b around the central frequency 
lol and to be off-resonant such that the interaction is 
well within the dispersive regime and absorption effects 
can be neglected. The detuning |A| is considered to be 
large compared to the Doppler width (^Doppicr and atomic 
decay rates Fatomic- Here and in the following, Fatomic 
denotes the largest effective atomic rate for transitions 
between ground state levels, including single particle as 
well as collective rates (see below). The magnetic field 
causes a Zeeman splitting of the atomic ground states 
fl. The strong y-polarized coherent beam is treated as 
classical field. With respect to quantization along x, it 
drives diagonal transitions |t) -^ |e^), \l) — >■ |e^). Figs.fl] 
and [2] depict only desired transitions, where photons are 
scattered into the copropagating x-polarized quantum 
field in two independent frequency bands, the upper and 
the lower sideband, centered around ujl ± ^• 

For the realization of the proposed scheme, several 
setups are possible. In a simple two-level model, where 
the Larmor splitting of excited states equals the Larmor 
splitting of ground states, a homogeneous static electric 
field can be applied to the second ensemble such that 
the resulting Stark shift enhances the energy difference 
between ground and excited states by 2A as shown in 
Fig. [2] This yields the following effective ground state 
Hamiltonian 



H = Ha 



Hr, + H\, 



where excited states have been eliminated under the con- 
dition |A| ^ Fatomic, fcoppier- Throughout the wholc pa- 
per, the convention h= 1 is used. Ha = ^ {Jx.i — Jxji) 
accounts for the Zeeman splitting of atoms in the exter- 
nal magnetic field and Hl = J dk (tuk — ojl) aj-flfc is the 



free Hamiltonian of the light field. In a rotating frame, 
the interaction Hamiltonian is given by 



Hint = / dk ^ g{k) U^ aj^.e' 
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-fdk y: .9(k) mE ^n.^^"' +-E -L^^^'^i«l 



+ H.C. 



(6) 



where the first and second integral cover narrow band- 
widths Aojis and Auj^s centered around the lower and 
upper sideband respectively. (A complete treatment 
based on the full Hamiltonian including all light modes 
can be found in App. |B 1[ ) Ak specifies the two orthogo- 
nal polarizations of the light mode with wave vector k. 
The atomic operator aj/jji — \t)i/ii,i{i\ refers to a par- 
ticle in ensemble //// at position r.^, Ak = k/, — k, and 
ki is the wavevector of the applied classical field. AC 
Stark shifts have been absorbed in the detuning. g(k)/i 
and 5(k)zy are the effective coupling strengths for the 
passive (beamsplitter-like) part of the interaction and 
the active (squeezing) component of the Hamiltonian 
respectively. More specifically, g(k)/i = -^z^dk and 

g()i.)iy = ^T'q" ffk ; where f^probc is the Rabi frequency 
of the applied classical field. Here and in the following, 
we refer to the off-resonant (entangling) light field as 
probe field. In the following sections resonant fields 
are introduced which will be referred to as pump fields. 
The definition of the coupling constant for transitions 
betw een g round and excited states g^ can be found in 
App 
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normalized such that /i 



1[57]. 



A Hamiltonian of type (|6| can be realized in many 
different ways. In general, the scheme presented here 
can be implemented in any system where a tunable 
quadratic interaction with an active and a passive part 
corresponding to two sideband modes can be realized, 
for example in ions or using optopmechanical resonators. 
We focus here on the creation of dissipatively driven 
entanglement in atomic ensembles and continue consid- 
ering the level structure depicted in Fig. [2] If the Larmor 
splitting of excited states is considerably larger than the 
splitting of ground states, it is for instance possible to 
introduce a A/2 plate between the two ensembles instead 
of applying an external electric field [S5j. As discussed 
in Sec. |IV[ alkali atoms provide another possibility to 
realize the desired light-matter interaction. Due to their 
multi-level structure is is not even necessary to introduce 
electric fields or passive optical elements. 

We remark for clarity, that the possibility illus- 
trated in Fig. [2] implies that the effective coupling 
constants (after adiabatic elimination) g(k)/i and g(k)z/ 
describing the interaction of light with the first and the 



second ensemble have different signs, as the light is blue 
detuned in the former and red detuned in the latter, 
such that /i/ = — /i// and vj = —vn- This is not the 
case in the implementation considered in Sec. |IV[ Due 
to the complex levels structure, both effective coupling 
constants have the same sign and therefore /i/ — fiji 
and ;// ~ vjj. In order to describe both alternatives in a 
compact way, we use a unified notation and absorb the 
sign in the definition of the atomic operators referring 
to the second ensembles un^i — >■ sgn(/i//i//)cr//.i as 
explained in App. |B2[ 



It is instructive to consider Hamiltonian (|6]), 
where excited levels have been adiabatically elim- 
inated, because it shows clearly that the light 
matter interaction depicted in Fig. [2] corre- 
sponds to a beamsplitting interaction of the type 

H « /a.,, ^k (Aai + Atak) + /^^^^ dk (Sai + B^a^) 
between photons in the upper and lower sideband with 
the nonlocal operators A and B (with additional phase 
factors e**'^'''"'). By deriving the corresponding master 
equation and including thermal motion as explained in 
Sec. |III C| one can show that this Hamiltonian yields a 
master equation which consists of a desired part of type 
(pi) with jump operators A and B and an additional 
contribution representing noise terms. However, in 
the following two subsections we derive the maser 
equation starting from the full Hamiltonian including 
excited levels since this approach is better suited to take 
dipole-dipole interactions into account. 



B. Effective master equation for ground states 

In the following, we outline the derivation of the 
master equation for atomic ground states |t) and \\.) 
and comment on the approximations used to obtain 
the shown result. The full calculation can be found in 
App. IT} 



For brevity, we use a short hand notation and 
abbreviate master equations of Lindblad form 
dtp{t) = k/2 {Ap{t)A^ - A'^Ap{t)) + H.C. with complex 
decay rate k and jump operator A by the expression 
dtp{t) = K/2Ap{t)A^ + ... . 

We consider the full Hamiltonian including excited 
levels and undesired transitions [59 without applying 
the rotating wave approximation for quantum fields (see 
Eq. Bl). As explained in App. Bl counter-rotating 



terms play an important role in the calculation the 
imaginary parts of the master equation, but do not 
affect the real parts. Starting from the full Hamiltonian, 
we obtain a master equation of Lindblad form for the 
reduced atomic density matrix (see Eqs. (B2)-(B4|). To 
this end, we apply the approximation of independent 
rates of variation [60| and follow the standard procedure 



assuming Born Markov dynamics. The approximation 
of independent rates of variations is valid if the Rabi 
frequency of the applied laser field flprobc is very small 
compared with the frequencies of atomic transitions. As 
we consider transitions in the optical domain, this as- 
sumption is clearly legitimate. More generally, here and 
in the following sections we consider situations exhibit- 
ing two very different time scales for variations in the 
system and in the bath of light modes FatomicTc << 1, 
where Tc is the correlation time in the reservoir. For 
optical frequencies this is very well justified and we can 
therefore assume Born-Markov dynamics. Moreover, we 
restrict ourselves to settings, where the level splitting fi 
between the states |t) and ||) is sufficiently large, such 
that the upper and lower sideband can be treated as 
independent baths [61 (compare App. [B1| ). Finally, we 
assume that the condition k^ ^ i?/L is fulfilled, fc^ 
is the wave vector of the applied laser field. Since we 
consider frequencies in the optical domain, k^ is on the 
order of lO^m"^. L is the spatial extent of the atomic 
ensembles, which we consider to be on the order of cm, 
while the distance between the two ensembles R is about 
one meter. 

As next step, excited states are adiabatically elimi- 
nated under the condition |A| ^ ratomici (^Doppicr- This 
leads to an effective master equation for atomic ground 
states. Using the abbreviated notation introduced in the 
beginning of this subsection. 



dtP{t)^\ Y. e-'^^'^^'^-^'h,, (A,p{t)A]+B,p{t)B: 
+ 2 E e-^'^^'-^'^J^i {c,p{t)CJ+D,pit)D 



«j'=i 



+ ... 



(7) 



where Jij and Jij are complex decay rates which are dis- 
cussed below and A = -j= X]i=i ^j- The operators B, C 
and D are analogously defined as sums. Ai, Bi, Ci and 
Di are given by 



A, 
B, 

D, 



= p, (7i,i + v a 



t 
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= P criiA + v cr 
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where the abbreviations cr^^j/jji = \ t)////,i(t I and 
<^UJ/ii,t = I Di/iiaH I are used. Terms involving the 
operators A and B represent desired transitions involv- 
ing a spin flip |t) -^ \i) or \l) -^ |t) as shown in Fig. [2] 
Terms involving the operators C or D represent unde- 
sired transitions which lead to dephasing [H]. Desired 
and undesired transitions are associated with different 
decay rates J^ and J^ respectively. In the four level 
model considered here, Jij = 2 Jij, due to the ratio of 



Clebsch Gordan coefficients c^ =±i/ca =o — 2. As in- 
troduced above, Jij = j{rij)+ig{rij) is a complex decay 
rate with real part 7(ry) — j{Yji) and imaginary part 
g(jij) = gijCji). Imaginary single particle terms represent 
energy shifts (single atom Lamb shift) and are absorbed 
in the definition of detunings. Therefore we consider in 
the following only imaginary terms gi^ij) with i ^ j and 
use renormalized atomic energies and the resulting effec- 
tive detunings. The real and imaginary part of Jij are 
given by [62] 



2\ sm{ki^rij) 

3„/, „,„ , .2\f cosjki^r^j) sin(fcLr^j) 
2'V- "- -'- )[ {k^r^^Y (fcLr,,)3 
3„/. ,^ . .2\ cos(fcLri-,) 



7(r,,) = -r(l-(p.%) 

|r(i-3(p-f,,) 



(9) 



2\i sin(fcL7-ij) cos(fcLr,y) 



(10) 



^T[i 3(p.f,,)^^ ^^^^^^^^ , ^^^^^^^^3 



where p is the unit vector of the dipole matrix element 
p — (e^|e x| t), which we assume to be real, f^ is the 
unit vector of the interatomic distance r^j = r^ — Vj and 
rij — Tji is the length of the vector Vij . T is the effective 
decay rate of a single isolated atom. 



C. Master equation including thermal motion and 
noise processes 



master equation is given by 



, N 

+ 2 Ef..(c.p(i)CJ-l-Ap(0^] 



ij = l 



+ 



(11) 



with Vij — 2Tij for the basic model discussed here. For 
fc,L»land^^j,r,,=rj^. d = iVJ- = j^ 

is the resonant optical depth of one atomic ensemble. 
Using this definition and 1/ (k^L) ^ 1, one obtains 



dtp{t) = d^Ap{t)A^ + d^Bp{t)B^ 



,r 



(12) 



4=1 

r ^ / 

+ ^^ 2 51 \^\,iP^^)^I,i + Al,iPi^)^II^ 



+ d^-Cp{t)C^ + d^-Dp{t)D^ + ^ ('^^ + ^') 

N 
'Y^{(^U,I,iP{t)(^U,I,^ + '7u,//,jp(i)cr;;,//,j) 



+ 



In this subsection, atomic motion is taken into 
account [551 - I55] . As is shown below, thermal motion 
gives rise to noise terms which are small compared the 
desired contributions for samples with high optical depth. 

Atoms are statistically distributed. The dynamics 
of the whole system is thereby governed by two different 
time scales, the characteristic time of radiative emission 
l/Fatomic and the characteristic time of atomic redis- 
tribution -, where L is the length of a cubic ensemble 
and V is the average velocity of particles. In the limit, 
where the time scale of atomic motion is fast compared 
to the time scale of radiative decay Fatomic^ "C 1 one 
can describe the emission independently of the evolution 
of atomic positions which enters the master equation in 
the form of averaged coefficients, where the average in 
time corresponds to an average in space |66j . Atomic po- 
sitions can be treated as independent random variables 
and for simplicity, we choose a Gaussian probability 

ESI. As 



distribution of width L, P{r) 



7r3/2L3' 



The first three lines correspond to the first sum in 
Eq. (11). The entangling terms in the first line are 
enhanced by a factor d. For sufficiently optically 
thick samples, additional noise terms in the second 
and third line, which reffect thermal motion, are 
small compared to the desired contributions. The last 
two lines correspond to the second sum in Eq. (Ill, 
where | t)(t I + I i)(l I = U was used. The first two 
terms d{r/2)Cp{t)C'< + d{t/2)Dp{t)D'< are collective 
dephasing terms. They do not have an effect on the 



shown in App. IBJ we find that imaginary parts of the 
averaged decay rates can be neglected. The averaged 



entanglement generated (see App. |C2[ ) and can therefore 
be omitted in the following. 

In the following sections, the effect of pump fields 
is considered. Resonant pump fields cause incoherent 
cooling (and heating) processes, which can be taken into 
account by adding cooling (and heating) terms which 
correspond to a transfer of atoms from level jj,) to level 
It) (and back). Finally, we include additonal processes, 
which do not lead to spin flips but cause dephasing, such 
as fluctuating magnetic fields. The full master equation 




FIG. 3: (Color online) Steady state entanglement ^oo in a 
two-level system versus Z = (|/i| — l^^l)"'^ for optical depth 
d = 30 per ensemble. The horizontal black line indicates the 
separable limit. (For separable states ^ > 1, the smaller ^, the 
higher the amount of entanglement) The lowest line (violet) 
depicts ^oo for purely radiative dephasing F^'*'' — 0. The next 
curves show in ascending order F^'*'* 
(green), F:^'*'* = lOF (red) and F^"^"* = 
is the single particle decay rate. 



= 2F (blue), Vf^ = 5F 
= 20F (orange), where F 



is given by 



(13) 



2 

J- heat 



N 

2 YlvL^P^^'>^i'^^'^n.tP{^)'^ii 



4=1 






Note that the last three lines represent single particle 
processes. They do not feature a collective enhancement 
factor d as the entangling terms in the first line. The 
noise terms proportional to F/i^ and Tv^ in the second 
and third line in expression ( |12[) have been absorbed in 
lines two and three of Eq. ( |13| . Hence, Fcooi (rhcat) 
is the total single-particle cooling (heating) rate. Noise 
terms proportional to F yp? + v"^) in expression ( 12 1 have 
been absorbed in the last line, such that Id is tlie total 
dephasing rate. More details concerning the derivation 



can be found in App. [C] For simplicity, we assume 
identical conditions for both ensembles. The amount 
of entanglement produced is measured by means of the 
quantity ^ defined in Eq. (l3|. Hence, the time evolution 
of Sj = var(Jyj -I- Jy.ii) + var(Jzj — Jz,ii) as well as 
the evolution of the mean value of the longitudinal 
spin \{Jx.i)\ — \{Jx,ii)\ need to be calculated. We con- 
sider the limit iV ^ 1 and start by computing the former. 

Tij decays according to 

dt^j = -(f + dTP2{t)\l:J + NdTP2{tf{\^i\-\v\f, 



The 



where f = F^ooi + Thcat + ^d and P2{t) = |(J^). 
evolution of the mean value of the longitudinal spin is 
given by 



dt{Jx 



' rt \ heat 



N 

rcool) \Jx}t + 7^ (rcool ^ Thcat) 



There are two distinct time scales. For atomic ensembles 
with high optical depth, the evolution of the transverse 
spin components is collectively enhanced and therefore 
fast compared to the decay of (Jx) which is due to sin- 
gle particle processes. In the limit where the entangled 
quantum state follows the changing atomic polarization 
adiabatically, the time evolution of ^(t) is given by 



m 



1 



-{r+drP2{t))t 



P2{t) 

1 f + dVB2{tY{\p\-\i^\f 



P2{t) f + dVP2{t) 

In the steady state 



(14) 
^ _ g-(f+drP2(«))A 



^o 



P-?. 



1 

-P2,oo 
1 cool 



dTPl^M 



iHr 



F 

^ heat 



dvp 



(15) 



2, DO 



cool 



heat 



This result shows that for high optical depth, the system 
reaches an entangled steady state. Under the dissipative 
dynamics considered here, entanglement persists for 
arbitraril y lo ng times. In the absence of noise, F = 
and Eq. ^[^ reduces to ^oo — {\p\ — 1^1) ■ 

Fig. [3] shows the attainable amount of entangle- 
ment in the steady state ^00 for moderate optical 
depth d — 2>Q versus Z = (|^| — |j^|)~ if only probe 
of the full master equation (fl3|) can be found in App. IC 2| fields are applied. In this case r^^oi'' = /^^E and 



^probe 

heat 



v^T. The dephasing rate Fd = F^ 



rad 



D. Creation of entanglement 

In this subsection, we determine how much entangle- 
ment can be generated by the proposed scheme in the 
presence of noise processes for a given optical depth d 
and given parameters /i and v. Details of the calculation 



2(m^ 



padd 
^ d 

i.2)F 



consists of a radiative part F™ '^'^° 
[57) . which is due to light-induced transitions |t) — > It) 
and II) — 7> II), and an additional term F^'^'^ which 
summarizes all non-radiative sources of dephasing 
such as fluctuating magnetic flelds. This additional 
component can take values up to F^'^'^ = 20F while 
still allowing for a reduction of ^00 by 15%. For large 
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FIG. 4: (Color online) Steady state entanglement ^J^* for 
optimal squeezing parameter Z versus pump parameter x. 
(Two-level model, d = 30). The curves correspond in as- 
cending order to Vf'^ = 15(grey), Vf"^ = 20r (orange), 
padd ^ 25r (brown), Vf"^ = 30r (pink) and Vf"^ = SST 
(light green). The inset shows the steady state polarization 
P2,oo versus Z in the absence of pump fields a; = (dashed 
line) and for a; = 5 (solid line). 



values of T^'^, the limiting mechanism is the decrease 
in polarization for high squeezing parameters and 
can be counteracted by applying resonant a+ and U- 
polarized pump fields to the first and second ensemble 
respectively, which drive the transition \\,) — >• je-t-). fn 
this case, the cooling rate can be roughly estimated 
as Fcooi = {l+x)TiJ? [68j. The pump parameter x is 



given by a; = 



')'lw 



(A-n)" 



k, where ilp 



is the Rabi 



-j=p— — n,, vviicic d!,pump 

probe 

frequency of the pump field and 7lw is the natural 
line width of excited levels. The correction factor k 
takes Doppler broadening due to thermal motion into 
account |69j . fn the presence of pump fields, radiative 
dephasing is enhanced V^'^ = 2 ((!-}- a;) ^^ + v'^) T. The 
heating rate is unaffected. Fig. |4] shows the maximal 
attainable amount of entanglement ^J^* (entanglement 
for optimal squeezing parameter Z), for d = 30 versus 
pump parameter x. For cc = 5, additional dephasing up 
to T'^'^ — 37F can be tolerated while still allowing for a 
reduction of ^J^* by 15%. Remarkably, the application 
of external pump fields, which amounts to adding extra 
noise to the system, is advantageous in this case. 



IV. 



IMPLEMENTATION IN MULTI-LEVEL 

SYSTEMS 



In Sec. Ill 



the theoretical framework for creating 
steady state entanglement between two atomic ensembles 
at room temperature is presented in detail for two-level 
systems. This section complements the main results 



derived in Sec 

in multi-level systems 



III by considering the implementation 
In the following we investigate 
the possibility of transferring the concepts developed 
for two-level systems to atoms with multi-level ground 



states by means of a general simplified model and 
analyze the conditions for obtaining entanglement in a 
steady state qualitatively. 

As specific example, we study the implementation 
of the proposed scheme in ensembles of alkali atoms, 
where the two-level system is encoded in a multi-level 
ground state manifold. Due to the richer internal 
structure, no external electric fields or optical elements 
need to be employed in contrast to the setup discussed 
in the previous section. As explained below, suitable 
values i-ii = flu and i^i — vn are realized naturally. In 
principle, it is possible to include all magnetic sublevels 
and all possible transitions of a particular alkali atom 
in the following consideration. However, rather than 
aiming for a complete description which takes the 
entire level structure of a specific atom into account, the 
general model used here is primarily intended to describe 
the underlaying physics. In Sec. |IV A[ we show how 
additional dynamics in a multi-level system can be taken 
into account by means of this simplified model which 
allows us to describe the physical effects with a small set 
of parameters while capturing all relevant features. In 



Sec. IV B we estimate the attainable entanglement. 



A. Including multi-level dynamics 

In the following we consider encoding of a two-level 
subsystem in a multi-level ground state. For example, 
the ground state of alkali atoms with nuclear spin / is 
split in two manifolds with total angular momentum 
F = / -h 1/2 and i^' = / - 1/2 respectively. The 
relevant two-level subsystem can be encoded in the 
two outermost states of the F = I + 1/2 ground state 
manifold |t) = \F,±F) and \i) = \F,±{F - 1)) in the 
first /second ensemble. 

In general, the maximum attainable amount of en- 
tanglement i^idoai _ (|/i| — |i/|) is determined by the 
different rates /i^F and i/^T at which probe-field induced 
transitions \l) — >■ |t) and |t) — > \i) occur. The values 
of the parameters /z and v depend on the multi-level 
structure of the excited states as well as on polarization 
and detuning of the applied laser field. This can be 
illustrated by considering the off-resonant probing of 
^^■^Cs atoms on the D2 line using the setup shown 
in Fig. [1] y-polarized probe light which propagates 
along z, interacts successively with two Cs ensembles in 
x-oriented magnetic fields. We assume weak magnetic 
fields {B w IGauss), such that the Larmor splitting 
fl ~ 300kHz is much smaller than the fine splitting 
of excited states. The first ensemble is strongly spin 
polarized along the orientation of the magnetic field, 
while the second ensemble is polarized antiparallel. As 
shown in Fig. [Sl the passive interaction, which transfers 
atoms from | |) to | f) involves the upper levels with 
F — 4, 5, whereas the active part of the light-matter 
interaction | f) ~^ \ V) involves the manifolds with 
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a) 




Ensemble I 



Ensemble II 



FIG. 5: (Color online) Implementation of the proposed 
scheme in ^'^''Cs ensembles. The atomic samples are placed 
in X polarized magnetic fields and interact with a y polar- 
ized probe field propagating along z as shown in Fig. [l] The 
relevant two-level subsystem is encoded in the 6S1/2 ground 
state, It) = |4, ±4) and ||) = |4, ±3) in the first/second en- 
semble and are coupled to the excited states in GPs/2- Only 
desired transitions are shown. Level splittings are not rep- 
resented true-to-scale. (In a magnetic field of 1 Gauss, the 
Zeeman shift of magnetic sublevels is about lO^Hz while the 
hyperfine, -and fine splitting is on the order of 10* Hz and 
10^* respectively.) 



F — 3,4,5. Taking the different Clebsch Gordan coef- 
ficients into account, one obtains Z = {p — i>) = 2.3 
for blue detuning A = 700MHz with respect to the 
F — 5 manifold of 6P3/2- (Both parameters, /i and v, 
are positive.) Alternatively, x polarized probe light can 
be used in combination with red detuning as shown in 
Fig. [6] In this case Z = 2.4 for A = -700MHz. Both 
variants are possible; in the following we will consider x 
polarized probe fields. 

More generally, the multi-level structure of excited 
states affects only the value of ^'deai_ rpj^^ multi- level 
character of the ground state leads to additional dy- 
namics that needs to be taken into account. Firstly, 
atoms can leave the relevant two-level subsystem. While 
matter and light interact, atoms redistribute or are lost 
to other groundstate manifolds. Therefore, the atomic 
population in the two-level subsystem is continually 
reduced. This is accounted for by introducing a time 
dependent population N2{t) and including the effect 
accordingly in the corresponding polarization P2(t). 
The subscript "2" emphasizes that these quantities are 

!>}■ 

we 
to 



defined with respect to the two- level subsystem { |t) , 



In order to calculate N2{t) and P2{t) in Sec. IV B 
introduce now a general model which allows one 
analyze the realization of the proposed scheme in 
atoms with multi-level ground states. A high degree of 
population and polarization with respect to the two-level 
subsystem is required in the process of generating 
long-lived entanglement. Therefore a± polarized pump 
and repump fields have to be applied. These additional 
fields induce transitions with Amp = +1 in the first 
ensemble and transitions with Amp = — 1 in the second 
one. For alkali ensembles, pump fields drive transitions 



within the manifold F ^ I + 1/2 while repump fields 
transfer atoms in i^' = / — 1/2 back to i^ = / -I- 1/2. In 
the desired case of high polarization with respect to the 
two outermost states, the atomic population in sublevels 
with F = A, ±mp < 3 in the first/second ensemble can 
be neglected. In this regime it is sufficient to restrict the 
description to three states, |t), \i) and \h) = \F',F') for 
the first and \h) = \F' ,—F') for the second ensemble, as 
shown in Fig. l6] 

Finally, one has to distinguish between spin opera- 
tors which refer to the relevant two-level subsystem and 
experimentally measurable quantities which are defined 
with respect to i^ = / + 1/2. For clarity, operators 
referring to the full multi-level structure are labelled 
with the subscript "exp". The longitudinal spin of each 
ensemble is given by 



N 



y y m,\m,)i{m,\ 



(16) 



i—1 m— — F 

N 



J2iF\t).{t\ + iF-l)\i),{i 



2F — 1 

J.,2 + -^^N2{t), 



(17) 



where N2 = J2z=i (It)i(t I + li>»U I)- For transverse spin 
components, 

N F 

J,,cxp = 2 E E VF{F + l)-mim+l){\m + l),(m| 

i—l 771 — -—F 

+ \m),{m + l\)^V2FJy^2, (18) 

such that 

Sj,exp = 2FSj,2 + 2(2i^ - 1)7V;, (19) 



B. Dissipatively driven entanglement between two 
alkali ensembles 

In the following, we outline the calculation of the 
entanglement which can be produced in the described 
setting and compute the time evolution of ^exp(i)- The 
master equation governing the evolution of the atomic 
system according to the general model outlined in 
Sec. |IV A[ as well as details of the calculation summed 
up below, can be found in App. iDl 



^N 



^N 



Atomic populations iV^ = Er=i lt)^(tl, ^;=Ei=i i;)^UI 
and Nil — J2i=i l^)i(^l can be calculated using the rate 
equations 



(20) 
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a) 




b) 


excited states\ 




lT>i 


^i 


^ \i) 




&3 \h) 



Ensemble I 



Ensemble 11 



FIG. 6: (Color online) General scheme for the realization of 
long-lived entanglement between alkali-ensembles. The po- 
larization of the probe field is assumed to be parallel to the 
applied magnetic field. The ground state S1/2 consists of two 
manifolds with total spin F = 1+1/2 and F' = I— 1/2, where 
/ is the nuclear spin. The states |t) and |4,) are encoded in the 
outermost levels of the F = I + 1/2 ground state manifold. 
We identify 1 1) = \F,F) and |i) = \F,F -1) for the first 
atomic ensemble and |t) = \F,—F) and ||) = \F,—F+1) 
for the second one. Only desired transitions are shown. In 
the presence of strong pump fields, the amount of entangle- 
ment generated can be estimated by a simplified three level 
model including the state \h} = \F' , ±F'), for the first/second 
ensemble. 



Sexp, 



0.96 



0.92 




X,. 



repump 



be- 



FIG. 7: (Color online) Steady state entanglement ^^ 

tween two ^'^^Cs ensembles versus strength of repump fields 



(see main text. Sec. IV B and App. ^ 



The plots 

depict data for ensembles with optical depth d = 30, blue 
detuned probe light with A = 700MHz and different values 
for the non-radiative dephasing rate Y"^'^ . The curves cor- 
respond in ascending order to T'"/'^ = (violet), T'f'^ = 2 
(blue), Vf'^ = 5 (green), Vf"^ = 10 (red) and Vf'^ = 20 
(orange) . 



where 



M 



(r-r; + Ttft.) 



r 



-t-i 
1ft 






ih) 






-2(r 



ir+r,. 



Tab is the single-particle rate for the transition \a) — > \b). 
If the transition rates are known, the number of atoms 
in the relevant two-level subsystem N2{t) = A'^ + iV^ 
and the polarization P2{t) = {N^ — N^) / N2{t) can be 



directly computed. Based on this result, the time evolu- 
tion of Sj,2 = var {Jyj + Jyji)^+va,T {J^j - J 2,11)2 '^^'^ 
be calculated. Again, the situation exhibits two different 
time scales. The decay collective of Ej^2 is fast compared 
to the evolution of N2{t) and P2{t)- A calculation anal- 
ogous to the one described in App. [C] leads to 



V,2 



Af2(0)e-(^+'*(*)r^^(*))* 



(21) 



N2{t) 



f + d{t)VP2{tf{tx-yY 



T + d{t)TP2{t) 

l_^-{r+d{t)TP2{t))t\ 



with d{t) = dN2{t)/N and f = % + T^^ + T^,, + T^, + 
Ttt+ri^+ri'^'^, where T^'^'^ accounts for non-radiative de- 
phasing. On time scales which are long compared to the 
fast desired dynamics (but short enough to avoid profuse 
depletion of the relevant two-level subsystem, such that 
N2{t) S> 1 is guaranteed) Ej.2 is given by the long-time 
(It) solution 



E% = N2{t) 



d{t)TP2{tf{^Ji~vY 

V + d{t)TP2{t) 



(22) 



Now, this result is related to experimentally measurable 
quantities. Inserting Eqs. (17 1, (18), (19) and (22) into 
the definition ^oxp = Sj,oxp/T2 I (J^oxpT) yields 



f" = 

Scxp 



f + d{t)TP2{tf{^i-v) 

f+d{t)TP2{t) 

N^{t) 2(2F-1) 



2F 



P2{t)+2F-l 



N2{t) P2(t) + 2F-l 



(23) 



This long-time solution is a generalization of Eq. ([5]) 
which takes multi-level dynamics into account. Particle 
losses result in a decrease in d{t) and P2{t)- If transi- 
tions out of the two-level subsystem can be counteracted 
efficiently by pump- and repump-fields, a true entangled 
steady state can be reached. Else, a quasi steady state 
is produced. These two cases are illustrated in Figs. [7] 
and [s] respectively using the concrete example of ^•^■^Cs 
ensembles (F=4) at room temperature. Fig. It] shows 
the amount of steady state entanglement generated in 
case of sufficient repump power. More specifically, the 
depicted curves represent solutions for different values 
of Fj'^'^ versus the repump parameter Xrcpump, starting 
from a^rcpump — 0.01. The repump parameter Xrepump 
quantifies the strength of the applied repump fields 
and is given by the ratio x^epump = ^?opump/^pump,opt' 
where f^pump opt is the optimal Rabi frequency that can 
be chosen for the pump field within the validity of the 
model considered here, {^p^^p opt is the minimal Rabi 
frequency of the pump field leading to N^,oo/N2,oo > 0.95 
in the steady state.) Details of the calculations leading 
to the plots can be found in App. |E] Fig. [8] iUustrates 
the case, where particle losses dominate the evolution of 
^oxp(0- The curves in this figure show the amount of 
entanglement generated in the absence of repump fields 
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as a function of time in ms. For short times the time 
evolution is governed by the desired dynamics within 
the relevant two-level subsystem and reaches quickly an 
entangled steady state. For longer times this stable state 
with respect to the entangling dynamics is superposed by 
the slow additional evolution imposed by the multi-level 
structure of the ground state. The fast desired dynamics 
entangles the collective spins of both ensembles, while 
particle losses cause a slow but continuing shortening of 
the spins. In this sense, Fig. [8] shows the creation of a 
quasi steady state. For t > O.OS/F, the time evolution is 
given by Eq. (23), that is, the time evolution is solely 



determined by particle loss-related dynamics. 

In this section, we put emphasis on the general, 
implementation-independent limitations of the proposed 
scheme imposed by radiative transitions. These unde- 
sired processes are characteristic for a given level scheme 
and intimately linked to the tradeoff between enhanced 
entangling dynamics due to increased probe, -or pump 
power and added noise. Depending on the concrete 
experimental realization, other undesired processes 
impairing the performance of the proposed scheme can 
occur like for example spin flips due to collisions. Atomic 
transitions and additional dephasing due to collisions 
with the walls can be taken into account by including 
terms of the type 



dtp{t) 



Fcol 



Fcol 
Icol 



1=1 

N 



i=l 

N 



} , {o'iiJAp{t)(Tiij^i + criiJiAP{t)'^UJL'' 



to the master equation. Since the thermal energy of 
atoms is typically much larger than the atomic level split- 
tings, one can assume the same collisional rate Fcoi for 
all atomic transitions. The value of Fcoi has to be deter- 
mined phenomenologically for the specific experimental 
setup under consideration (compare [56]). 



V. CONCLUSIONS 

In conclusion, we propose a technique for entangling 
two mesocopic atomic ensembles at room temperature 
which are separated by a macroscopic distance. The core 
idea is to engineer the coupling of the atomic system to 
its environment in such a way that the steady state of 
the dissipative time evolution is the desired inseparable 
state. As entanglement is produced in the steady state 
of the system, it is long-lived and immune to noise. 
The reservoir consists of the common modes of the 
electromagnetic field and the coupling of the bath to 



Cexp(0 




FIG. 8; (Color online) Quasi steady state entanglement of 
^'^''Cs ensembles versus time in units l/F in the absence of 
repump fields (a;rcpump = 0; all other parameters take values 
as in Fig. ml. The fast entanging dynamics results in a drop in 
Ccxp(i) (indicating the creation of an inseparable state), but 
since particle losses are not counteracted by repump-fields, 
this additional slow dynamics eventually causes Cexp(i) to rise. 
Hence, a quasi steady state is produced. The behavior on 
long time scales is determined by the stationary state given 
by Eq. ( 15 I superposed by slow multi-level dynamics. 



the system can be controlled by means of laser- and 
magnetic fields. We provide a detailed theoretical 
description including dipole-dipole interactions for two- 
level systems and find that the imaginary parts of the 
master equation (collective Lamb shifts) are negligible. 
Hence, light-induced collisions do not play an important 
role in the setup considered here. The proposed scheme 
for generation of entanglement by dissipation is analyzed 
for two level systems and complemented by considering 
the implementation in multi-level ground states. 

Future directions include the transfer of the ideas 
presented here to other systems, where a quadratic 
interaction involving two sideband modes can be real- 
ized. The system may either be described by bosonic 
modes (compare Sec. O and App. pi) or spin degrees 
of freedom. If the Hamiltonian corresponding to the 
interaction between this system and light can be de- 
composed into an active and a passive part such that 
each part involves one sideband, as in Hamiltonian (|6|, 
dissipatively driven entanglement can be generated using 
the method described in this paper. 



Optomechanical resonators [8TH83j interacting with 
light are promising candidates. If an optomechanical 
system is driven by a strong pump laser, the linearized 
radiation-pressure Hamiltonian gives rise to a passive 
(beam-splitter) interaction for positive detuning be- 
tween cavity and pump- frequency [84j [85] . The resulting 
effective interaction is active (squeezing) for negative 
detuning. The effective optomechanical coupling rates 
can be adjusted by tuning the intensities of the driving 
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fields [55]. Hence, the quadratic interaction in this 
system provides naturally the basic prerequisites for 
implementation of the proposed scheme. One can 
for example envision a setup, where light interacts 
subsequently with two movable mirrors which each 
constitute one end of an optical cavity such that the 
nanomechanical oscillators are driven into an entangled 
state. Equivalently, membranes coupled to Fabri-Perot 
cavity modes could be used. 
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Appendix A: Steady state entanglement for bosonic 
modes 



In the following, it is shown that the two mode 
squeezed state ptms = |*tms>(*tms|, with A|*tms) = 
^I^TMs) = 0, is the unique steady state of the time 
evolution described by the Master equation 

dtp{t) = KA (Ap{t)A^ - A^Ap{t)/2 - p{t)A^A/2\ 

+ Kb (Bp{t)B^ - B^Bp{t)/2 - p{t)B^B/2\ . 

As stated in Sec. [Hj the bosonic mode operators a and 
h {[a,a^] = \b,h^] — 1) can be transformed into the non- 
local operators A and B by the unitary operation 

A = UaU^ = pa + v b\ 
B = UbU^ = pb + iy aK 

Since unitary transformations preserve commutation re- 
lations, [A,A^ = [-8,-8'^] = 1 with ^2-1/2 = 1. By insert- 
ing these expressions in the equation above and defining 
pu — W pU we find 

dtpu(t)^KA{apu{t)a^ -a^apu{t)/2-pu{t)a'^a/2) 

+ KB{bpu{tp-b^bpu{t)/2-pu{tpb/2).iAl) 

This is a master equation for two modes coupled to a 
bath with temperature T = 0. The steady state is the 
vacuum |0,0)(0,0|, witha|0,0) =6|0,0) = [Zg. Hence, 
inverting the unitary transformation yields the unique 
steady state C/|0,0) = |*tms)- 



For bosonic modes, the amount of entanglement 
can be quantified by means of the violation of a local 
uncertainty relation in terms of quadratures [TTJ [72] . 
For entangled states 

var(a;+) -f var(p_) < 1, 

where x+ = (xa + Xh)/V2, P- = (pa - Ph)/V2 and 
Xa = {a + a^) /^/2 and Pa — —i{a — a^^ /\/2 (analo- 
gous expressions hold for x\y and pb)- bi particular, 
var(a;+) -I- var(p_) = [p — v) = e^^'', for two mode 
squeezed states with squeezing parameter r. 

For large, strongly polarized atomic ensembles, col- 
lective spins can be described by bosonic modes 
7^E»=\cr/,^ « a, ^7=ESo-//,* ~ ^ using the 
Holstein-Primakoff-approximation |73j . In this case, 
^ < 1 (see Sec. Ill| is equivalent to the criterion 
var(a:;+) 4- var(p_) < 1. 



Appendix B: Derivation of the master equation 

In this appendix, the master equation for creating 
entanglement between two atomic ensembles at room 
temperature discussed in Sec. |III[ is derived in de- 



tail. App. |B 1| and App. |B2| complement Sec. IIIB and 
Sec. |III C| respectively. In the former, we derive the mas- 
ter equation for atomic ground states, Eq. (B5|. In the 



latter, we include thermal motion of atoms and calcu- 
late the effective decay rates. We show that in the setup 
under consideration, the resulting master equation can 
be assumed to be real, since imaginary parts play only a 
minor role. 



L. Master equation for atomic ground state levels 

It) and l;) 

Light and matter are assumed to interact as described 



in Sec. HI A We consider the full Hamiltonian including 
undesired transitions 59J and without applying the rotat- 
ing wave approximation for quantum fields. It is given by 
H = Hi, + Ha + Hint, where Hi, = J dk {ujk - ujl) al-Ofe 
is the Hamiltonian of the free light field and if a accounts 
for atomic energies in the rotating frame. H^ — Ha.i + 
Ha,ii with Ha,i = J2ii^tj\(^t)i:ii^t\ + ^lj\(^l)iA'^l\) 
and Ha,ii = J2zi^tJi\(^t)ii,i{(^t\ + ^lJI\(^l)II,i{(^l\)■ 
}leTe, we introduced the detunings A^ //// and Ai^j/jj 
which correspond to diagonal transitions \\.) — >■ |e-|-) and 
It) ^ |e^) respectively in the first/second ensemble. In 
the setup illustrated in Fig. [21 A^j = — A^j/ = A — fi 
and A^ / = — A^ // = A -I- IL The interaction Hamilto- 
nian Hint — Hci + i?qu consists of a classical part Hd, 
which accounts for transitions induced by the driving 
field and a quantum part -ffqu, which involves quantized 
field operators. The Hamiltonian describing the interac- 
tion of light with the first atomic ensemble is governed 
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by the Hamiltonian 



(Bl) 



N 



H,i,i = r!p,obe V e'^-'^' {h)iAi\ + H)iAt\) + H.C., 



i=l 



N 



N 2 



3ke''"''«kae;)z,,a|e'^" + |e^),,,(t|e-'"* 



+ |t)/,.(e^|e'%-2»<.^*)+iJ.C., 



,5ke''"-'ak(|et)/,.a|+|e;)z..(t| 



+EEE. 

where Ak specifies the two orthogonal polarizations of 



«=1 k xt,=l 



the light mode with wave vector k. .gk = ^k ■ Pi / 5717 

is the coupling strength of desired transitions involving 
the quantum field, while gk is the coupling strength 
corresponding to undesired transitions, ek is the unit 
polarization vector, V is the quantization volume of the 
electromagnetic field, eg the vacuum permittivity and p 
is the transition matrix element of transitions |e^) — > | l) 
and |e;) -J' I t)- -^int = gint.i + gint.ii- Hmtji is given by 
an expression analogous to Eq. (Bl|, where subscripts 
are changed accordingly (/ -^ II). 



Based on this Hamiltonian, we derive a master 
equation for the reduced atomic density matrix pit) us- 
ing the approximation of independent rates of variation 
and assuming Born-Markov dynamics as explained in 
SeciniB] 

dtp{t) = -i[IIcup[t)] + CquP{t), 

where £ is a Lindblad operator corresponding to the 
quantum part of the interaction Hqn- This master equa- 
tion consists of three parts 



dtp(i)=L(p(t))ens.I + i(/O(t))cns.II+i(p(0) 



intcr-cns. 



.(B2) 



The first and second term L(p{t))cns.i and i(p(i))cns.ii 
include only terms referring to the first and second en- 
semble respectively, while the third term i(/5(i))intcr-ens. 
summarizes all terms combining operators acting on both 
samples. The first term is given by 

L{pit)),ns.i = -i[H,ij,p{t)] (B3) 

+ ^E'^<"lt)/.^(^tlpW|eT)/,,(tl 

i,3 



where we used the approximation w^ ^ fi, which is very 
well justified for optical frequencies, and neglected fast 



oscillating terms which appear in the standard derivation 
[60] , if photonic modes in the upper and lower sideband 
are not treated as independent baths [ST] . J^- = J-^ is 
a complex decay rate associated with desired transitions, 
is the rate of undesired transitions. 



while J^^ 






III 



rl-J 



For the simple model discussed in Sec 

Imaginary single particle terms Im J^^ ) (Lamb 

shifts) are absorbed in the detunings as explained in 
Sec. IIIB The second term L(p(t))cns.ii is given by 
an analogous expression with changed subscripts / 1— )■ //. 



The last term in Eq. (|B2|) can be written as sum 

Hp{'t)YintLcns.- The 



^(P(^))inter-ens 

first part is given by 



'^(P(*))intcr-cns. 



Lip{t)) 



I. II 



^E^''''lt>/.(eTlpW|e;>//,,ai 

i,3 

^E^''''lt)/.(e;|pW|eT)//,,ai 

id 

(B4) 



«j 



\ii,i 



and ^(p(i))int'or-cns ^^ givcu by an analogous expression 
with changed subscripts (/ — > //) . Decay rates appearing 
in inter-ensemble terms differ from single-ensemble rates 



J. 



I.I 



II. II 



J. 



n-IJI _ -ikR^IJ 



^11, 1 _ +ikR ^ 



1,1 



where R is the distance between the two ensembles. In 
order to obtain compact expressions, we use the simpli- 
fied notation J'ij for single ensemble or inter-ensemble 
rates depending to which samples the indices i and j re- 
fer. We use the convention Yij = Tj — r^ — R if the atom 
with index j is located in the second ensemble. (How- 
ever, in App. |B 2[ it is shown that the distance between 
the ensembles does not play a role in the setting consid- 
ered here.) Using this notation, J'ij is given by 



J,, =ydk^(f (k)e'''(--'-'-^)y d-ic 



g i(t^i;-WL)T I g-i(wfc+wi,)-n 



Ak = l 



where the sum over light modes was changed into an 
integral J2k ~^ T^^/c^k- The prefactor is absorbed 

in the coupling constant g(k) = ^JVji^lTxY ffk whenever 
an integral over light modes is used. The second term 
in the expression in brackets (e~*("'=+'^i)'^) stems from 
counter-rotating terms in the Hamiltonian, and would 
not appear if the rotating wave approximation had been 
applied. Using the identity l^ e"^"^ = 7r5(a;) -I- 17^(1/0;), 
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where V is the principal value, we obtain 

2 



Re{J,,)=7rJdkY^g\ky 



■ik(r^— Tj) 



(5(wfc- Wi), 



Ak = l 



Im{J,-)=tvijdkJ29\^)e'''^'''' 



Ak = l 



1 



1 



lul - Wfe ujl + i^k/ 



These rates can be calculated as shown in [S5] (compare 
Eqs. ([9]) and ([10|). 



Now, a master equation for the reduced density matrix 
of atomic ground states Pgp(i)Pg is derived by apply- 
ing the projector P^ = (g)^^! (|t)7,«(tl + li)/.*(il) ® 
(lt)//,i(tl + l4-)//,j"(4-l) to the differential equation 

dtp{t) = L{p(t))ens.I + i(p(i))ons.II + -^(p(i))intcr-cns. 



using Eqs. (B3) and (B4|. Excited states are eliminated 



under the condition A^.////, A^j/// >> ratomic, using 
dtFep{tWg = dtPgp{t)P, = dtFeP{t)Ve = 0, whcre 
Pe = U — ¥g. Moreover, we assume that terms corre- 
sponding to states with two or more excitations, for 
example terms of the type Pe 'p{t)Pg, are negligible 
compared to terms corresponding to states where at 
most one atom is in an excited state like Pi p{t)Fg. 



with 



Pe'^ = EtlPe,/,. (E> Vgji + E^^.Vgj (g> Pe,//J 

Ve,^J/II = (8)ti(|et>////,.(eTl + |ei)////,.(eil). Pi'^ 
is defined analogously. In the following we denote the 
resulting reduced density matrix of atomic ground states 
Fgp{t)Pg simply by p{t). We obtain 



dtpit) 



n 



probe y^g-iki,(rj- 




where the abbreviations (J^^j/n^i — \ t)////,i(t I ^-^^d 
<^U,i/ii,i ^ I i)i/nA{i I "^ere used. Terms with pref- 
actors 1/A'^ have been neglected since we assume that 
detunings are large. AC-Stark shifts 



N 



dtp{t)\ 



AC Stark 



Zi iprobc / ^ 



^Z7 



a 



H;I ,i 



^\,I 



Pit) 



N 

^^ 'probe / ^ 
1=1 



zf2r 



'^n,II,i I '^U,II,i 



Hj 



A 



t.i 



, Pit) 



are absorbed in the detunings. Using the definitions 



A-i-n 



A-n 



571 and Ji. 



^^/" - *^7a!!' ^^/" - -2VAn 
Jij 2r2pi.obeVAr2/ (A^ — Q'^), one obtains 

dtp{t) = - J2 e-^^^^'-'^'^J^j (Ap{t)A] + B,p{t)B}) 



N 



1 J2 ^^-^^-^^U,, (c.p{t)C] + D,p{t)D 



i,j=i 



+ ... 



(B5) 



with 



A^ 


= Pi Oui + Vu cr}^ ,, 


B, 


= Pll OlIA + Vi <t\^^, 


a 


= P-I f^UJ-i + '^11 f^ttJIyi 


D^ 


= PlI f^U-II-i + ^-f '^ttJ,i 



The expressions in the main text are obtained by intro- 
ducing a unified notation (compare Eq. ([8])), which ab- 
sorbs the relative sign pi/pn in atomic operators refer- 
ring to the second ensemble an — )• sgn(/i//i//)cr//. 



2. Master equation including atomic motion 

Below, we include thermal motion of particles, by 
treating atomic positions as classical random variables. 
We start from the master equation for atomic ground 
states (B5|. As outlined in Sec. IIIC random atomic 



positions can be taken into account by introducing aver- 
aged coefficients in the master equation. Averaged rates 
( Jy ) are calculated assuming Gaussian distributions with 
width L for atomic positions in the two ensembles. We 
start by calculating the rate corresponding to moving 
particles in a single ensemble. Below it is shown that for 
inter-ensemble rates the same result is obtained for the 
setup and range of parameters considered here. 

iJ'^j) = ^ /dr/dfe*'^-(--^)-^(7(r - r)+zg{r - r)). 



7r3L6 



(27r)^/^ L3 



dr_e 



iit-Lr- 



^(7(r_) + ig(r_)),(B6) 



where the variable transformation r+ = r + f , r_ = r — f 
was made in the second step. The single particle rate 
is given by Ja = T (Lamb shifts are absorbed in the 
detunings). Now, averaged rates (Jij) with i ^ j need 
to be determined. We consider now first the real part 
and then the imaginary part of the averaged decay rate 

The real part Fy- — RejJJij)) is calculated by in- 
serting Eq. ^ into Eq. (B6) and fixing p = x and 



^ij,A + ^ij.B is a sum of two contributions 
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corresponding to the first and the second line in Eq. Q. 
The first term is given by 



3r 



ij.A 



dr^r_i 



"■- sin(fcLr_) 



\J'1-kL'^ Jo k^r- 

I sin(fcLr_) cos(fcir_) sin(fcir_) 



(B7) 



fcrr_ 



{kLT-Y 



[kLT^r 



The integrand of this expression tends to zero in the 
Hniit r_ — > 0. This can for example be seen by 
expanding the integrand for small values r_ <C 1, 

^ sin(/cir_) /sin(/cir_) |^ cos(/i:i,r_) sin(/i;i,r_) 



r_e SL-i 



2^2 
3 - 



^(, 



cos(A;i,r_ ) 



0{rt_). Hence the dominant contribution in the 
limit kLL 3> 1, stems from the first term in brackets 
sa\{kLr-)/(kLf'-). The other two terms in brackets de- 
cay faster in the interatomic distance r_ and lead only 
to corrections on the order of l/(fciL)*. For fc^L ^ 1, 



3r 



- ij,A 



1 



' 9 -^ sin(fcLr_) 



{kLT-Y 



4 (kLL) 



, (1 - e-'^-^) 



which can be approximated by TijA ~ |r , .a . The 
second part 



ij,B = 



3r f°° 2 -^ cos(kLr) sm (kir)' 
ar-r_e 21,2 ' 



(kLrJ) (kLT-) 

3cos{kLr} 3sin(fcLr_)\ 




(kLr-Y 



(kLT-Y 



is negligible compared to the first part ^ij^A in the 
asymptotic limit fc^L — )• 00. Its integrand tends 
to zero for r_ — > (expansion for r_ ^ 1 yields 
A;|ri/45 + 0{r^)) and contains only terms proportional 
to cos{kLr-)am{kLr_)/ (kLT-Y, cos(kLr-Y /(kL^-Y 
and sm(kLr-Y /(kLf-Y with x > \. These types of 



terms have been neglected in Eq. (B7| or decay even 
faster in r_. Since Yij,B is negligible compared to r^.^, 
weuser,, =r^^^. 

Next, we calculate the imaginary part Gij = Im({Jij)) 
by inserting Eq. (10) into Eq. (B6). As before, we 
consider the two contributions Gij = Gij^A + G ij,B 
corresponding to the first and the second line in Eq. ( 10 ) 
separately. The integrand of the first part 



G 



3r 



I'}, A 



\/2^L3 JO 

sin(fcLr_) 

kLT- 



dr-r_e 2l 

cos(/cLr_ 
(kLT-Y 



- COs(kLT-) 



kLT- 
sin(/cir_] 
(kLT-f 



tends to zero for r_ — > (expansion for r_ ^ 
1 yields 2r_/(3fci) + 0(r^Y ^^d features a rapidly 



oscillating term proportional to sin (fc^r-) cos (fciT-) 
in the integral, which leads to a contribution which 
scales with l/(/ciL)''. The other terms proportional 
to cos^ {kLL)/(kLLY and cos{kLL)sm{kLL)/(kLLY are 
again of the type discussed and neglected before. Hence 
it is well justified to assume that Gijj ^ ^ij.i- The 
integrand of the second part 



Gy,B - 



3r 



27ri3, 



dr-r_e 2t2 



sinfe^r-) cos (fci,r_) 



(kLT-) 



(kLT-Y 



s\n.(kLr-) 3cos(kLT-) 3sin(fcir_) 



kLT- 



(kLT-Y 



(kLT-Y 



also tends to zero for r_ — > (expansion for r_ ^ 1 
yields —r-/(15kL) + 0(r'^)) and contains only one term 
which has not been considered so far. The term in the 
integrand proportional sin(fcir_)^/(fcir_) leads to to 
a contribution which decays with \og(kLL)/(kLLY in 
the asymptotic limit fc^L — >■ 00. The imaginary part 
of the averaged decay rate (Jij) is therefore negligible 
compared to the real part. 

The distance between the two atomic samples does 
not play a role in the setting under consideration. In the 
limit fci ^ R/L^, averaged single ensemble rates equal 
averaged inter-ensemble rates (Jij) = {Jij}^'^ = {Jij}^'^^- 
In the following, we outline the calculation of the inter- 
ensemble value F- : ^. Analogous arguments can be used 
compute T^-g, G^^^, and G^^g. 

Inter-ensemble rates are obtained by averaging atomic 
positions with respect to Gaussian distributions centered 
at the origin and a distance R apart respectively 

^ fdr fdfe'^'^^''^^'>el^ e~^'L-^^\(r - f), 



p7,//_ 



7r3L6 

-X^fdv+ldv-e'^^'' 
ttYl^J J 

1 



(2 



-(._+R,)2 -(r+-R)2 

e 21:2 e 2Z2 7(r_), 



3/2 



(2ttY' L^ 



dr-e 



ikiF- 



-(r_+R)2 



7(,r-J 



where the variable transformation r_|_ = r + f , r_ 



r — r 



was made, as before. By inserting Eq. ^ and neglect- 
ing the dipole factor (x • (r — f)/|r — f |) , which does not 
play a role for the distance under consideration |74| , one 
obtains 



^i.ii 



3F 



2(27r)^/^L3 
3F 



, »k^r -(--t^)' sin(fcLr-) 

dv- e ^ e 21:^ — ; , 

kjr_ 



""dr-r'jMfil^rdeMe) 



2V2ttL^ Jo "' ' kLT- jg 

gjfci,r_cos(e)p^(r^+-R,^+2r_flcos(e)) 

3F 

"" /o 



2V27rL3fcL 



6 21,2 

dr- sin (kLT-) e 21^ 



ei"' 
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In the last step, the integral over 6 was transformed using 
the variable transformation x = — cos(0)r_. The integral 
over X can be directly evaluated yielding 



r 



/,// ^ _ 3r 



/•CO 

/ dr_ sin (/clT-) 

"'0 



e 2L^ 



e 

3r 



ik^v- 



-ikLr--\-- 






/oo 
dr_ sin {k^r^ 
-oo 



e 2t 



^ — g«fct'-_ 



As next step, the variable transformation f = r_ + i? is 
made such that 



r 



I, II 



3r 



2V27rL3/cL (ifci - i?/L2 



3r 



AV2^L3kL {-kL - iR/L^) 

ar 



(if sin {kL{f — i?)) , 



dfe 21-2 ^ 



AiL^kl + ikLR) 



1 - e~"'=^ 



2k; L^~2i 



kLR\ 



which yields r[j^^ = ^iTfc^IF ^°'^ ^^ ^ ^/^^ ^^"^ kLL^ 



Appendix C: Time evolution of entanglement in a 
two-level model 



In this appendix, we calculate the amount of entan- 
glement produced ^(i) (compare Eq. (Is])) for the model 
described in Sec. |III[ The first part of this appendix |C 1| 
covers the derivation of ^{t) based on the full master 
equation (13). The second part C2 contains explana- 



tions concerning Eq. (131 



1. Time evolution of entanglement 

In the following, the time evolution of entangle- 
ment ^(i) is calculated. To this end we calculate the 
single-ensemble variance of transverse spin components 
var(J2) — {JD — {Jz)"^, the inter-ensemble product of 
transverse spins (JzjJzji) and finally the mean value 
of the longitudinal spin (Jx)- 

We start by calculating dt{J^)t- The dissipative 
evolution described by Eq. ([l3| leads to 



Applying the decorrelation approximation (JzJx) ~ 
{Jz){Jx) [Z5] for mean values of products of transverse 



and longitudinal spin yields 
and similarly 



N /- 



r + drP2W> +^^ 



N /- 



dt{Jl) = - (f+drP2(i)) {jD + ^ (f + dTP2{t)\^Ji^ + V 

where (J^) = f P2(i) and (Jj) « {J^f = ^P2{tf were 
used. The latter approximation leads only to an error 
of the order O (;^) [75] Next, the mean values of the 
transverse spin components are computed using the same 
approximations. 



dt{Jy/z)t = -]^{f + dTP2{t)){Jy/z)u 

with {Jy/z)t=o = 0- The mean values can therefore 
be ignored when calculating single ensemble variances 

var(Jy/^). 

The time derivatives of inter-ensemble products of 
transverse spins are given by 

dt{Jz,IJzJI)=-{T + dTP2{t)YJzJJz.II) + ^^iydTP2{tf, 

dt{JyjJyJl)=-(r + dVP^it)) {JyjJyJl)-^fll^dTP2itf, 



where we used (JxjJx.ii) ~ {Jx,i){Jx,ii)- For N ^ 
1, this is a very good approximation, since collec- 
tive effects on populations are suppressed by a fac- 
tor d/N. Therefore, the time evolution of longitudi- 
nal spins is only determined by single-particle terms, 
which do not lead to correlations between the two en- 
sembles. Hence, the variances of the non-local operators 

Jy,± = {Jy,I ± Jy.Il) / V^ and Jz^± = {Jz,I ± Jz,Il) /V^ 

evolve according to 



dtvar(Jy^±) 



f + drPaW var(Jj,^±) (CI) 



^(f + drP2(t)2(MT^f 



dtvai{Jz,±) = -(f + drP2(i) var(J,,±) (C2) 



^(f + drP2(t)2(;i±i.f), 



such that 

var(Jy^±)^ = var(J^^zp) 



NT + dTPl^itiTi^f 
4 f + dTP.^ 



in the steady state. The variances var(Jj,.+) and 
var(Jj,__) are squeezed, while var(Jj,_) and var(J2_|_) 
are anti-squeezed. Now, we consider the time evolution 
of the longitudinal spin 



dt\{Jx) 



N 



(Fheat — rcool) — (Fheat + rcool) l(>-'a:)|, 
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which yields directly 



I , ^ , I _ N Fcool — Thcat _ ^ p 

I Wo;/ loo — ~;r t; ~t; — ":7-f2,co 

^ ^ cool "r i heat ^ 



for t — >■ c». 

Collective effects have an negligible effect on the 
time evolution of the polarization. P2it) evolves 
due to single-particle effects only and hence much 
slower than var(Jj,-|-) and var(J2^±) for samples with 
high optical depth. In this case, the solution for 
^(t) = var {Jz,+) /P2{t) can be cast in a simple analytical 
form 



m = 



1 



P2{t) 
1 f 



-(f+drP2(t))t 

dTP-2{t)\\il\ 



P2{t) 



T + dTP2{t) 



_ -(r+drP2it))t 



2. Full master equation 



In the following, we comment on the form of Eq. ( 13 1, 



in particular on the absence of collective noise terms. 



The probe fields considered in Sec. Ill are off-resonant 
and it has been shown in the main text that collective 
contributions feature an enhancement factor which ren- 
ders them the dominant decay mechanism for samples 
with high optical depth. As is shown in Sec. HID 
it can be advantageous to apply also resonant laser 
light (pump fields). In contrast to off-resonant fields, 
collective contributions are negligible compared to 
single-particle terms for resonant light in the situation 
considered here. Unlike off-resonant collective rates, 
resonant collective rates are much slower than the 
corresponding single-particle rates for samples with 
high optical depth, which is an effect well known and 
harnessed in electromagnetically induced transparency 
[77H5D] . The single particle decay rate after adiabatic 

elimination of excited states is given by Tres — ''"'"'' , 
where fipump is the Rabi frequency of the applied laser 
field and 7lw is the natural line width of excited levels. 
Coherent effective effects lead to an enhancement factor 
d in the denominator. Intuitively, this effect can be 
understood by noting that emitted resonant photons are 
reabsorbed in an optically thick medium. 

Due to atomic motion in ensembles at room tem- 
perature, spectral lines are Doppler broadened. We take 
therefore off-resonant contributions of pump fields to 
the master equation with a detuning on the order of 
the Doppler width fcoppicr into account. A calculation 
along the lines of the derivation shown in Sec 



CI 



shows that these terms are negligible compared to their 
single-atom counterparts. More specifically, collective 
terms corresponding to a detuning feoppior lead to decay 



rates proportional to """■"'' 7lw'^ = ^^^^^ f -n^wd \ 
with |(5Doppicr| >> 7Lwd, while single particle resonant 
terms lead to decay rates proportional to — esse. 



Finally, Eq. ([13|) 
corresponding to 



does not include collective terms, 
radiative processes which do not 



change the internal atomic state, since they do not have 
an effect on the amount of entanglement generated. A 
master equation corresponding to the terms omitted in 
Eq. pi 



dtp{t) = d^Cp{t)C^ + d^Dp{t)D^ + ... , 

with operators C = J2i (m f^UJ,^ + '^ ^n.H-i) ^"^^ ^ = 
^ = J2i{^^ ^UJi.i + ^ ^nJ.i) leads to dtCit) = 0. Since 
MJy) = -§j{Jv), dt{J.) - -|^(J.), and dt{J,) = 0, 
(Jy) ^ (J^) = and (J^) = N/2 for all times. The 
time derivatives of single-ensemble variances for trans- 
verse spin components is given by 



dtij'y) 

dt{Jl) 






such that {Jy) — {JD ~ -/V/4, for all times. Accordingly, 
(Jy.iJy.ii) = {JzjJz.ii) = for all times since 



dt{JyjJy.ii) = -2—Tp,u{JzjJzji), 
dt{JzjJz..ii) = 2 — tp,v{JyjJyji). 

The processes under consideration do not create entan- 
glement unlike the terms in Eq. (12) with jump opera- 
tors A and B. As shown above, they do not degrade 
entanglement either. Collective terms corresponding to 
far off-resonant radiative transitions | t) — >■ |e^) — )■ | t)j 
I 4-) ^ |e^) — > I i) do not introduce random phases and 
preserve coherence. The emitted photon does not reveal 
information about the internal atomic state, since it is 
emitted into the laser mode. Terms with jump opera- 
tors C and D lead only to very small correction terms 
proportional to 1/iV and can be ignored. 



Appendix D: Generation of steady state 
entanglement in alkali atoms 

In this appendix, we consider the generation of 
dissipatively driven entanglement in multi-level systems 
based on the model described in Sec. IIV Al 

Taking three ground state levels | t), | |) and \h) 
into account, as shown in Fig. [6] the evolution of the 
reduced atomic density matrix can be described by the 
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master equation 

dtp{t) = dTApit)A'f + dTBp{t)B'< (Dl) 

N 



i=l 

N 



i=l 

N 

rthE(l'^)^^'(t|pWlt)/,^(/i|+|/i)//,.(t|pWlt)//,^(/^l) 
1=1 

AT 

T/^E i\t)iAh\p{t)\h)i,,(^\ + \DiiAh\p{tmiiA\) 

i=\ 
N 

r^E ih)iAi\pit)\i)iAh\ + h)HAMm)iiAh\) 

i=\ 

N 

T^Y. i\i)iAh\p{t)\h)iA\ + \i)iiAh\p{t)\h)iiM 

i=\ 

N 

TtrE in)iAt\p{m)iAt\ + n)iiAMm)iiAt\) 

i=l 

N 

r^^E i\i)iAMm)iAi\ + \i)iiAMm)iiAl\) , 



i=l 



where Tab is the single-particle rate for the transition 
|a) — > |6). As in Sec. |IIIB we omit collective terms due 
to resonant pump fields, as well as collective dephasing 
terms (see App. C2). Collective terms involving the 



level \h) are also insignificant for A^ >> 1, as long as the 
number of coherent collective excitations is small. 

In order to compute the amount of entanglement 
produced, we consider the variance of the nonlocal op- 
erator Jj/,+,2 = {Jy,i + Jy,ii)2 /A^- (The subscript "2" 
emphasizes that these quantities are defined with respect 
to the two- level subsystem {It);!^)}-) A calculation 
analogous to the two- level derivation in App. |C 1| shows 
that {Jyj)2 = {Jy,ii)2 = for all times. Therefore 
var ( Jj,)2 = {Jy)2- For simplicity, we assume that both 



ensembles are identical {Jyj)2 = {Jy.ii)^ 



{Jl)2- 



According to Eq. (|Dl|), the time derivative of the 

72\ 



single- ensemble variance {Jy)2 is given by 



dt{ry)2 = - r 



dT^m) 



iJt 



.N2{t) 



where the decay rate F, the number of atoms in the 
relevant two- level subsystem (It)?!-!')}! A^2(0 and the 



corresponding polarization P2(t) are defined in Sec. IV B 



N ^J2^ (jt)»(t I + \i)t{i I + \h)^{h\) is the total number 
of atoms in one ensemble and A^2(0) = N. Note that 
repump fields, which transfer atoms from \h) to jf) or 



I4-) (cor resp onding to terms with prefactors T^g and Ths 
in Eq. (Dl |) do not contribute to F. 



Inter-ensemble correlations {JyjJyji)2 evolve according 
to 



dt{Jy,iJy,ii)2 = - F + dr 



dV 



.N2it) 



N 

1^2(i)„,,,2, 



P2{t)]{Jy,lJyJl)2 



2 N 



-P2{tYpV. 



Hence, the time evolution of y&v {Jy^A)^ 
{Jy,iJyji)2 is given by 



dtUl- 



Analogously, 



dtiJ'h 



dT^P2it)){JlA2 



dV 



N 

1^2(i)„,,,2 



4 N 



P2{tY{p-^Y 



{J'y)2 + 

4 
(D2) 



N 
c^F-^^P2(t) (p-i^) 



P2it)]{JlA2 + r 



N2(t) 

4 
(D3) 



Since the evolution of N2{t) and P2{t) is known from 
equations (20), S]j,2 — {JyA^ + (Jz-)^ ^^^ ^"^ directly 
calculated yielding a complicated expression. However, 
as explained in Sec. HI and Sec. IV N2{t) and P2{t) can 



be considered to change slowly compared to the fast en- 
tangling dynamics. In this case, (D3| leads to the simple 



and convenient expression (21) used in Sec. IV B 



Appendix E: Implementation in hot 133Cs vapors 

In this appendix, we apply the results derived in 
Sec. |IV| to a specific example and consider the generation 
of entanglement between two ^•^■^Cs ensembles at room 
temperature. The parameters used in the following 
take values consistent with the experiments reported in 
P^ [50] . The approximate calculation outlined below 
provides a rough estimate of the entanglement that can 
be produced. 

We assume y-polarized probe light which propagates 
along z and interacts in succession with two ensembles 
in a magnetic field which is oriented along x. The laser 
field is assumed to be blue detuned by A = 700MHz with 
respect to the 6Si/2{F = 4) — )• 6P3/2{F = 5) transition 
(D2 line). Fig. [5] depicts the relevant parts of the atomic 
level schemes in both samples and illustrates the atomic 
transitions due to the light-matter interaction induced 
by the applied laser field. Initially, all atoms are pumped 
to state It). The restriction of the analysis to the three 
levels It), li) and \h) in the presence of strong pump 
fields, as described in Sec. IV A for x-polarized probe 
light is also valid for this configuration, as the rates 
of transitions from level | f) to states with iriF = ±2 



19 



occur at rates which are two orders of magnitude 
smaller than transitions within the sub-system under 
consideration. (r|4^4)_j.|4_2) 

r|4,4>^|3,2) = 0.02 r|4 4^^|4 3^) 



= 0.03 r|44^_^|4_3^ and 



In order to calculate the experimentally measur- 
able steady state entanglement using Eq. ( 23 ) for a given 



optical depth d and parameters /i and v, -^2(^)7 P2{t) 
and r need to be computed. N2{t) = N^{t) + A^;(t) 
and P2{t) = (iV^(t) - Ni{t)) /N^it) are readily obtained 
from Eq. pO| ), if the rates for all transitions are known. 
As the probe field is assumed to be off-resonant, probe 
induced rates Fab for transitions \a) — >■ \b) are calculated 

using the formula T^™ "" = flp^obo YjI a^^^lw, where 
the sum runs over all levels contributing to a particular 
transitions (for example the states |5,3), |4, 3) and 
|3,3) in e^Pg/s if T^; = Tv'^ is computed), /^f is the 
detuning for each contributing level, 7lw the natural 
line width of excited levels and cf° is the product of 
the corresponding Clebsch Gordan coefficients. Pump, 
or repump induced transitions are resonant and involve 
to a good approximation only one level. Hence the 



corresponding rates are given by T 



pump 
ab 



tii 



and r™P"™P = "lw""* '^^cpump k respectively, where 
k = ^Doppicr = aim ED- Here, we consider pumping 
resonant with respect to the level |4, 4) in QP1/2 (Di line) 
and pump fields resonant with respect to the level |4, 4) in 
6^3/2 (D2 line). (The decay rates are approximately the 
same ^Oi ~ ID2 = 7lw-) Having expressions for N2{t), 
P2{t) as well as f = % + Tj.^ + Tt?, + T^a + % + % + r 
at hand, the amount of entanglement which can be 
produced in this particular setting can be calculated. 
Results are shown in Figs. [7] and |8] and are discussed 
in Sec. IV B As the three-level description becomes 
inaccurate if too many atoms are transferred from 
state 1 1) to state 1 1) , both plots show results for the 
optimal (that is minimal) pump power which guarantees 
a fraction of at least 95% of all atoms in the relevant 
two-level subsystem in state |t) for all times. Besides 
the need for sufficient pump fields, repumping of atoms 
from P = 3 to -F = 4 is required. If strong pump- but 
no repump fields are applied, no entangled state can be 
reached, as for t — >■ cx), all atoms are transferred to level 
\h). 



7LW pump 
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